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Abstract 

In a very recent paper[J. Chem. Phys. 128, 114308 (2008)], Springall and co-workers reported 
Quantum Monte Carlo calculations on the potential energy curve of the helium dimer. They 
argued that their FN-DMC results are mostly within 1% of the Hurly-Mehl curve[J. Res. Natl. 
Inst. Stand. Technol. 112, 75 (2007)]. In this comment we show that their results are wrong. 
Their reference-results computed from Hurly and Mehl's curve are completely inaccurate, and also 
their FN-DMC results are wrong and unauthentic. 
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In a very recent paper Springall and co-workers reported Quantum Monte Carlo cal- 
culations on the potential energy curve of the helium dimer. They employed Slater-Jastrow 
form trial wave functions and used the fixed node approximation for the fermion nodal sur- 
face. In their work, they present the results of variational Monte Carlo VMC and fixed node 
diffusion Monte Carlo Carlo (FN-DMC) calculations for the helium dimer with atomic sep- 
arations in the range R= [0.9,7.4] a.u. Their FN-DMC result for helium-helium interaction 
energy at R=5.6 a.u is very accurate, predicted -10.89±0.17K or -10.96±0.15K, which are in 
very good agreement with the current accepted values of around 11.00K. But, it seems very 
strange that their other results at different separated distances are significantly different 
from the well published results. For separated distance at R=4.5 a.u, the well published 
and accepted values are all around 58.40K(Selected results are listed in TABLE 
but their FN-DMC result is 49.76±0.29K, there is about 9K error; for the separated dis- 
tance at R=4.3 a.u, the corresponding well published and accepted results are all around 
118.0Kpj-|ls|, but their FN-DMC result is 102.24±0.32K, this is nearly 16K error compare 
to the well published and accepted results. These errors are very large in such accurate 
Quantum Monte Carlo calculations, the errors are compare to, even large than ground date 
interaction energy around 11.00K at R=5.6 a.u. It is impossible and unbelievable that a 
benchmark method(Quantum Monte Carlo) could give a accurate interaction energy at one 
distance(R=5.6 a.u) but give inaccurate results(so large errors) at other neighbor distances. 

To investigate the possible error in their calculations, we carefully performed systemic 



Quantum Monte Carlo calculations for helium dimer [20] . As is shown in the TABLE II, our 
results are in excellent agreement with other well published and accepted results. Even more 
surprising, in their origin paper, their listed results computed from Hurly and Mehl's 21 1 



analytic potential for the helium dimer are also significantly different 



results. We then carefully re-computed the Hurly and Mehl's [21] analytic potential at 



rom the well published 



several distances. We find that Hurly and Mehl's 



2l| analytic potential should lead to 



two $ slightly different potentials, as they are shown in TABLE II. These two potentials 
are significantly different from the author's results, but, they are very good consistent with 
the most recently SAPT based analytic pair potentialjs] for the helium dimer and other 
theoretical results {4, 11 |. So, the results computed from Hurly and Mehl's 21] analytic 
potential in their paper are inaccurate. And near all of their FN-DMC results (except 
R=5.6 a.u) are very different from these analytic potentials and well published theoretical 
results. All the aboving results are listed in TABLE II. 



In summary, we conclude with pointing out main error in Springall and co-workers's 
original paper. They incorrectly computed analytic potential for the helium dimer from 
Hurly and Mehl2l|. The fixed node Diffusion Quantum Monte Carlo calculations in their 
paper should be wrong. Thus, their results are in very error, it was proposed to avoid, their 
Quantum Monte Carlo simulations should be revised, a more careful and responsible study 
is needed. 

Acknowledgement: We thank Prof. Krzysztof Szalewicz for their fortran program for 
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TABLE I: Comparison of selected predictions of the helium-helium interaction at R=4.5 bohr. 



Energies in Kelvin. 



Method 


Interaction energy 


Fixed node DMC a 


49.76±0.29 


Green function QMC[3] 


60.0 


Analytical potential[4] 


60.44 


rl2-MR-ACPF [5] 


58.49 


CCSD(T)+FCI correction[6] 


59.54 


SAPT [7] 


58.037 


CCSD(T)[8J 


59.470 


SAPT [8] 


58.371 


CCSDT(Q)/CBS[9] 


58.397 


CCSD(T)-R12 [10] 


59.543 


ECG [n] 


58.517 


CCSD(T)/CBS+FCI [12] 


58.407 


EQMC [13] 


58.3 



"The Springall and co-workers's FN-DMC result 
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TABLE II: Interaction energy values of the helium dimer at different values of the atom separation. 
The interaction energies are given in K. 



R(bohr) FN-DMC a Hurly b Hurly c Hurly d SAPT FIT e FN-RMC^ Other results 



0.9 


314409.15±3.79 


319211.08 


348590.96 


348588.32 


350152.17 


350482.26 


1.9 


45959.52±4.11 


45840.26 


44863.43 


44862.40 


44847.24 


44852.74 


2.3 


20159.33±0.64 


20108.76 


18728.15 


18727.51 


18716.36 


18718.39 


2.8 


5575.28±5.37 


5549.73 


6021.75 


6021.39 


6020.67 


6021.23 


6.1 


zz I l.y4±1.4o 


zzd9.1o 


1666. 6 1 


qqqq i r: 

zooo.15 


qqqq ro 

zooo.58 


QQQ A 1 1 

zoo4.il 


3.8 


537.78±1.77 


535.64 


508.66 


508.55 


508.42 


508.75 


4.3 


102.24±0.32 


100.67 


118.13 


118.07 


117.93 


118.01 


4.5 


49.76±0.29 


50.538 


58.547 


58.502 


58.406 


58.47 


4.7 


19.49±0.27 


21.020 


24.312 


24.276 


24.220 


24.229 


5.6 


-10.89±0.17 


-11.05 


-10.9957 


-11.0092 


-11.0048 


-11.003 


5.9 


-10.31±0.15 


-10.40 


-10.1793 


-10.1893 


-10.1865 


-10.184 


7.4 


-2.57±0.18 


-3.412 


-3.3297 


-3.3322 


-3.3330 


-3.332 



350489.6075 
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results, to be submitted . (see: http://arxiv.org/abs/1001.3268l 
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